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6. The angle increases io° a year, 

7. Mag. from Cor, D.M. Estimated 80. Is 87 in C.P.D. and marked: 

Angle decreasing, distance decreasing. 

There is a 13th mag, star 200° ±, 35" ±. 

8. Change. Two careful estimations with the chief star out of the field. 


Short Method for the Calculation of the Orbits of Celestial Bodies . 

By D. A. Pio. 

(Communicated by the Secretaries .) 

1. Gain of Time is the Object. —The determination of the 
elements of the orbits described by a celestial body has been 
carried by Gauss to the highest pitch of perfection the present 
state of mathematical science allows us to reach. The author 
has not the pretension to propose a new method for solving the 
same problem ; he only intends to substitute short, simple, and 
easy calculations for those which the method of Gauss requires, 
and limits himself to the most simple case : how to determine 
the first approximate values of the elements of a newly dis¬ 
covered celestial body with the least loss of time, so as to render 
the calculation of a new orbit an easy performance. 

On the authority of the illustrious French mathematician 
M. H. Poincare, the writer begs to state that the calculation of 
the first approximate orbit, by Gauss’s method, requires from 
fifteen to seventeen hours. A reduction to about six to s$ven 
hours is what the author is working for. A gain of about eleven 
hours with a sufficiently approximate determination of the 
elements Would be a service to practical astronomers. 

The author proposes, for attaining this short calculation of 
orbits, a modification of Gauss’s method containing innovations 
of great practical interest for the efficient simplification of the 
mathematical work. 

The writer begs to remark that the problem of orbits has 
been really solved long ago; that the question is only how to 
render tedious calculations more easy, but not to present new 
theories by which the work of La Place, Olbers, and Gauss may 
be swept away. This is the point of view from which the author 
begs his work to be considered. Practical improvement, not 
radical subvertment, is the only amelioration Gauss has left 
posterity to accomplish, 

2. Preliminary Remarks. —It is supposed that the data fur¬ 
nished by the observations and those taken from the solar tables 
are corrected for precession, nutation, and aberration of ohe 
fixed stars. Parallax and aberration in time are not taken into 
consideration. The calculations are made with logarithms to 
not more than five decimal places. 


© Royal Astronomical Society • Provided by the NASA Astrophysics Data System 


Downloaded from http://mnras.oxfordjournals.org/ by guest on August 25, 2015 





aKi* ■ t?9' •syhnw£06i 


Dec. 1903. Orbits of Celestial Bodies, 135 

3. Simplifications of Gauss’s Method. — (1) Gauss established 
the equation 

3 ft . sin 4 z = sin (z -f w) 

which leads to the value of r 2 . This part of the solution cannot 
• be improved. However, modern authors have brought confusion 
into the signs (+ or —) of the many auxiliary quantities 
occurring in the calculation : especially the case where r 2 is 
less than unity offers difficulties. The writer has reverted to the 
, Theoria Motus , and has removed all ambiguities (see articles 4 
and 5). 

(2) The solution of the equation 

3 ft. sin* 2 =; sin (z + w) 

has been rendered shorter by using logarithms with three deci¬ 
mals for the first approximations. The method of logarithmic 
differences has been reduced to its simplest form (see articles 7 
and 8). 

< (3) Gauss deduces the elements of the orbit from the known 

values of r x and r 3 , and from the angle 2f = v 3 ^v z contained 
between them. The calculation of these quantities requires the 
previous calculation of the distances from the Earth, p T and p 39 
which as done in the usual manner is the most lengthy and the 
most tedious part of the whole work. The usual formulae, as 
given in Tisserand, Legons sur la Determination des Orbites, 
pp. 76 and 77, are twenty-four in number, and some are very 
troublesome. The writer has reduced them to six by using the 
“ fundamental equations ” for the calculation of x x =p x cos J 3 Z and 
r 3 = p 3 cos fi 3 . The time thus gained amounts to some hours 
(see article 9). 

(4) The parameter, the eccentricity, and the longitude of 
, the perihelium of an orbit are the most difficult points in the 
' determination of the elements. Gauss invented a most ingenious 
method for the calculation of the parameter by successive approxi¬ 
mations through auxiliary tables. The author has much shortened 
the calculation of the parameter by reverting to the formula of 
Euler given in Theoria Motus , § 86 (VII.)* It is the shortest 
existing. 

(5) The formulae usually employed for the calculation of 
the angle of eccentricity, and the true anomalies v z and v 3 are 
very tedious. They occasion a considerable loss of time (see 
Tisserand, p. 80). The arc v z has been determined by a simple 
formula taken from Theoria Motus . The value of the eccentricity 
may be calculated very easily when the parameter, one of the 
radii vectores, and the corresponding true anomaly are known. 
The author, profiting by this circumstance, calculates two sepa¬ 
rate values for e and takes their geometrical mean as the most 
probable value. 

Each of the above modifications is certainly a detail of no 
theoretical importance by itself, and has no other merit than 
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136 i¥r. Pio, Calculation of the 

either to shorten the calculation or to increase the accuracy. 
The author begs to submit that all these modifications, taken as 
a whole,* give quite a new appearance to the method of Gauss. 

The question is, then, whether the above modifications satisfy 
or not the requirements of practical astronomers. 

4. Rules of Signs .—The calculations preliminary to the solu¬ 
tion of 

9ft . sin* z = sin (z -f w) 


are relative to the determination of the quantities S, r, T, 2, 
rr , S, d , u , w, and h. The angles ^ r, and u are to be taken 
always as positive arcs less than 18o°. The arc <j is always less than 
90° : it may be positive or negative, but it must bC the least of all 
values that result from the absolute number found for tg (<r + 1 ) ; 
and the sign of tg (2 - £ + rr) must be determined in conformity. 
The arc w must always be less than 90°, but it may be negative 
Vor positive. The signs of the linear quantities T, S, d , O, and h 
are -f- or — according as the calculation determines them. 
Every other rule, even given by great authorities, leads to 
impossible results. 

Note. —Oppolzer’s rule to make T always positive is false 
when r 2 is less than 1. In this last case w also is positive and 
not negative. 

5. A Difficulty in the Solution of 9ft . sin* z = sin (2-f-«).— 
The equation 

9ft. sin* z = sin (z -f w) 


is usually applied only to the planetoids where r 2 is comprised 
between 1*5 and 4, and then w is usually negative and nearly 
equal to z. In such cases the unknown value of z is found by 
successive approximations, the value of — w being the first, and 
the method of logarithmic differences leads very rapidly to a 
yalue sufficiently approximate for the first determination of the 
orbit. However, when r 2 is less than H 2 , then 8 is greater than 
90°, and the angle w is not only positive but very different from 
z, so that it becomes difficult to find a first approximation for z. 
As usually w is small, by assuming sin u == o the equation 


becomes 


and then 



9ft sin* w = sin (z-f it>) 


9ft. sin3 z = cos to 


sins: = 



cos <«> 



which may serve as first approximation. » ' 

6. Logarithms to three Decimal Places .—In the usual manner 
calculating the orbits much time is wasted by using logarithms 
%> many decimal places, while the Results foupd out ‘with so 
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much labour are scarcely trustworthy to the third figure. This 
happens with the first attempts to solve the equation 

2$. sin* % = sin (z + w) 

Also the correction of the approximate value of the parameter 
(Euler’s formula) needs only logarithms to three decimal places. 

7. Logarithmic Differences .—The method of logarithmic dif¬ 
ferences given by Gauss in Theoria Motus , § 11, for the solution 
of Kepler’s Problem is excellent for solving the equation 

2ft , sin* % =: sin (z + w) 

Its shortest form is that the author proposes. 

Example.—^[ 1*0822321] , sin*£ = sin(s—12°), where 2 ft 
= N[i ’0822321], w = —12 0 , A = logarithmic difference of sin z 2 
for one minute (respectively one second), p = log dif. for i f (or 1”) 


of sin (z 2 + (x)) } A = log sin* s 2 — 

log sin* z t , z t = 

- U). 

log sin z x 9’3 l8 

4 

X = o-5 

log sin* z 2 9*362 

4 

log sin* z x 7*272 

fi = $-6 

log sin* z 9 7*448 

log 1*082 /i 

1 

II 

Cn 

A 

-- 49 ' 

log sin (« 2 + hi) 8-354 

A = 176 


**+<<> = i° if'5 

0 

M 

M „ 

1 

II 

3 

«* = i 3 °i 7'-5 

*3 = r 4 ° 6'-5 


Six lines with logarithms are sufficient. 

8. Calculation of p x . cos / 3 1 and p 3 . cos f 3 3 . —The condition that 
the centre of the Sun must lie in the plane determined by the 
three observed positions of the celestial body furnishes three 
“ fundamental equations,” (5), (6), (7), in Tisserand, p. 47, con¬ 
taining the unknown quantities = p t . cos / 3 „ r 2 = p 2 . cos / 3 2 , 
r 3 = p 3 . cos / 3 3 . One of them, r 2 , being known, the determination 
of ti and r 3 becomes very easy. Equation (7) furnishes directly 
Ti , and from (6), where r x and r 2 are known, r 3 is easily found. 

9. Concluding Remarks. —(1) The abridgments introduced 
by the author will furnish the values of the elements with suf¬ 
ficient accuracy, when the intervals of time between two successive 
observations are not more than ten days for planetoids. 

(2) The purpose of a first determination of an orbit being 
first the calculation of a provisory ephemeris, and next the 
setting-up of numerical values for the elements approximate 
enough to allow the easy correction of these numerical values by 
comparison with a series of observations, any luxury of calcula¬ 
tions, as those indicated in the text-books, is a mere waste of 
time. Whatever be the care bestowed upon the first determina¬ 
tion of an orbit, the result arrived at by so much work will be 
discarded . Therefore brevity ensuring sufficient accuracy is the 
ideal of the first determination of an orbit. 
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(3) The calculation of an orbit is considered as the trial- 
piece of an incipient astronomer* It is, indeed, a proof of pro¬ 
ficiency, especially if carried out by the method of Gauss, as 
explained in Oppolzer’s text-book. But even the dreaded 
method of Gauss may be popularised, and that is what the 
author has attempted. In about six to seven hours, when the 
observations have been already reduced and corrected, to deter¬ 
mine with sufficient accuracy the numerical values of the 
elements of an orbit without employing special auxiliary tables, 
only using the common logarithmic tables of numbers and 
trigonometrical functions to no more than five places of decimals, 
this is certainly no achievement in theory, but only a sort of 
mechanical craft. However, is it useful or not ? Has the writer 
succeeded in building up an engine of calculation proper for the 
use of practical astronomers ? \ T 

10. A Test for Asteroids .—The writer has chosen as an 
example for the application of the simplified method of Gauss 
the, same observations from which Gauss himself deduced, the 
erbit of Juno in Theoria Motus , § 150 to § 155. His calcula¬ 
tions are a model of accuracy, and a comparison with those made 
by the writer shows the degree of precision the first values for 
the elements may reach when the modifications, proposed by thd 
author are used. The following table contains in its first column 
the results of the “ last hypothesis ” of Gauss : 



Gauss. 

Pio. 

Differences. 

r 2 

2*1183 

2*1187 

+ 0*0004 

2 f 

7 ° 34 ' 54 " 

7 ° 35 ' 7 " 

+ I $" 

a 

' 2*6450 

2*6499 

+ 0*0049' 

e 

°‘ 2 453 2 . 

0-24744 

+ 0*00212 

i 

13° 6' 44" 

13 0 7 ' 3 6 " 

+ 46" 

a 

.171° 7' 49" 

171 0 8' 29" 

• + 4 P" 

n 

K 52 0 18' 9" 

5 2 ° 30' 23" 

+ 12' 14" 

m " 

824"*8o 

00 

ts> 

tii 

O 

- 2"-3 

X 

41 0 52' 22” 

4t° 58' 22 " 

' + 6 f o" 


Formula. 

1. Symbols .—For the three observations of the planet at the 
epochs E t , E 2 , and E 3 , we have : 

Ej = t T . 

E^-— E 2 === t 2 . ) 

R I} R 3 , R 3 , radii vectores of Earth. 5 

L„ L 2 , L 3 , geocentric longitudes of Sun. 

* A x , A 2 , A,* geocentric longitudes of planet, 

d n 'fi„ fa, /> 3 , geocentric latitudes of planet. 

/ : : ^/ 2 , ^ heliocentric longitudes of planet. 

b J9 b 2) 63, heliocientric latitudes of planet. 1 
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r„ r 2 , r 3 , radii vectores of planet. 
v 19 ' v 2 , v 3 , true anomalies of planet. 

j0 2 , p 3 , distances of planet from Earth. 
e, eccentricity of planet. 
p, half parameter of planet. 
a , half major axis of planet. 

\ T d , planet’s revolution in mean solar days. 

m", its mean diurnal motion in seconds of arc. 
i, inclination of planet’s.orbit. 

S3, longitude of planet’s ascending node. 

IT, longitude of planet’s perihelium. 

L 0 , mean longitude of planet at epoch E 0 . 

E 0 , epoch of departure. . 

2. Preliminary Calculations .— 

A = R,[ty/ 3 3 . sin (X,—Lj)—ty/ 3 ,. sin (\ 3 —L,)], (i) 

B = R 3 [^p 3 . sin (A,—L 3 )—ty/ 3 , . sin(A 3 -L 3 )] ... (2) 

‘ C = R 2 [^/ 3 3 :‘sin(A I —L 2 )—'<$ 3 ,. sin (X 3 —L 2 )] ... (3) 

D = —U,[tgft 3 . cos (\ 1 —'L?)—k/S 3 I . cos (X 3 —L 2 )] ... (4) 

to.\p = -— / q ^ 2 ‘ positive and <180° ... ..." (5) 

sin(X 2 —L 2 ) 

tgc = ± ^( X °~ L 0 . (6) 

. cos-'/' ... 

Note.—< >90° when (A 2 —L 2 )^ ?°o j 2 < 9 °° when (A 2 —L 2 ) 

> 90 p. 

<270°’ 


2 positive and < i8o°. 

tgr = —--- r positive and < 18o° 

y . sin (A 3 —Aj). 1 



sin t 


to 2 =. ±... 

y T.sin(r + ^) 

Note. — a always <Cgo°, positive or negative, 
or negative, so that a has always its least value : 


•• - X7> 

.. ... ( 8 ) 

... (9) 

... (10) 

Take 2 positive 


S = 


C. sin 2 


sin 


P = 




d= 


P+i 


Ad-PB 


S 


... (u) 

... (12) 

- 03) 
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C?. COS . cr 
&7<7 . COS 


. Mr. Pio, Calculation of the 
u positive and <180° ... 


comp, log V 2 = 9*84949 - 


tg<D — - _ * 

v 2 . sin (45°—%) 

Note .—id is always <90°, positive or negative. 

^ __ d . sin <r 
sin w 

h= 1}(A + PB)R 2 3 sin3 S 

f, = J *,(*! +2 Q . 

U — Kt*(2tt + Q log \ — 5-69301-10 
o o 


m 


_Aj z + FBj 2 
“ h 


3. Solution of 2D£. sin* z = sin (z -f w). 

4. Calculation of r x and r 3 .— 

R 2 sin ?> 


r Q = 


sin % 


__ R 2 . cos / 3 2 . sin ( 5 — z) 


sxn« 


Y _ «? x 

'“r 2 3 - 

y 2 = h... 


r* 3 


= 




n. 


t x +t 2 

_ 


(i+YO 


~(<+Y 2 ) 


+ 4 

r x = r 2 . sin (X 2 —X 3 ) + w x R x sin (L x - X 3 ) t 

+ w 2 R 3 sin (L 3 —X 3 ) — R 2 sin (L 2 —X 3 ) f 

: n x sin (X x - 

r r Jg/ 3 2 -n 1 x 1 tgl 3 1 - - 
3 3 


ctyc x = 


r x 


R x sin (X x —L x ) 

:= ^1 + £1 


cot (X x — L x ) 


lxiv. 2, 
... (14) 

10 (15) 

... (16) 

... (17) 

... (18) 

... (19) 

... (20) 

... ( 2 !) 
... ( 22 ) 

— ( 23 ) 

... (24) 

••• ( 25 ) 

... (26) 

A 3 > ( 27 ) 

... (28) 

... ( 29 ) 

- ( 30 ) 
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“‘< X ’ «. 

... (30 

.l 3 = X 3 .+ ? 3 ... , ... ..... 

... (32) 

tgl sin^-Xd . ‘ 

y 1 R x sin (\ x —L x ) 

- (33) 

, fl ;X-_’ : 3 ^ 3 8in -( / 3 - X 3 -) . 

* 3 E 3 sin (X3-L3) . . 

- (34) 

*. _** fix 

r i . , ••• ••• 

sin o x 

- (35) 

r 

3 sin b 3 

... (36) 

5. Calculation of Inclination and Ascending Nods- 


,<</&!. sin 

tgb ? —tgbi . cos (/ 3 —h) 

... (37) 

D D< i 8 o° or > i 8 o° 

Note whether the motion is direct or retrograde. 

... (38) 

tgi = . 

smD 

-• .(39) 

For direct motion £ <90° ; for retrograde i > 90°. 


6. Calculation of Eccentricity and Major Axis., — 


• / v sin (L ~ If) . cos . cos b 3 

Sm (tU ^—V---- 4 ... 

.. cos^ , 

... (40) 

\/ p ' = r '%Xl) lo g^ — ?' 2 3 5 5 8 ... 

... (41) 

... 

... (42) 

^ = cot (v 3 —v z ) — rAp—rj)^ - _ 

^(p-^sin (v 3 -v,) 

,.. ( 43 ) 

e ,_ (P~rx)(p-r 3 ) 
r-sX 3 cos i>i. cos 1/3 

... (44) 

<*= p 2 ... ... ... ... ../ 

1 —e 2 

( 45 ) 


T* = a* X 36s-2s6 d log 365-256 = 2-56260 (46) 


• , *$" = log 1296 000= 6-11260 (47) 
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LXIV.J2, 


7. Calculation of the Longitude of Perihelium and of Mean 
Longitude .—- 

__ tg~D 


tgu 1 


COSfc 


(48) 


When D <180° Wi <180° ; when D >180° u x > 180°. 


n —— + 8 -“Vi mi <4* •••. (49) 



M = w —e.sinw; ... ... ... ... ... (51) 


M and w to be taken in parts .of the radius. 

r‘ i =E I —E 0 ... .(52) 

L" = + M x 206 265" 

.log 206 265 =5*31442 ..(53) 

Lo = L+n ... ••• #•* ••• (54) 


Calculation of Juno 9 8 Orbit 

1. The Data .— ‘ 

4 = 11^96324 4 = 9 d *97ii9 

log = 9*99968 log R 2 = 9*99810 log R 3 = 9*99697 
L x = 192 0 28' 28" L 2 = 204° 19' 49" L 3 = 214 0 16' 10" 

= 354 0 44' 32" X 2 = 352 0 34 ; 22" X^ = 351° 34' 

ft = -4° 59' 3 i" ft = -6° 21' 55 " ft = —7 0 17' 51" 

v . >.. / 

2. Preparatory Calculations. — 

f . log P = 0*07910 ^+^ = 21^93443 

log 4 = 1*07785 log 4 = 0*99875 log(4 + 4)= 1*34113 

log (4 + 24) = i* 5 ° 3 8 7 l°g (24-M2) — i* 53 OI 7 

Xj—L t = 162° 16' 4" Xi —L 2 = 150° 24' 43" 

Xj — L 3 = 140° 28' 22 /r 

X 3 — L,; = 159 0 6' 2 V x 3 — La = 147 0 14' 41" 

X 3 — L 3 = 137 0 18' 20" 

X 3 -X, = — 3 0 10' 2" X3-X, = - 59' 52" 

X 2 — X x = — 2 0 10' 10" X 2 —L 2 = 148° 14' 33" 

Lg—L a = 9° 56' 21" L 3 —L x = 21° 47' 42" 

L a —L, = n 6 51' 21" 
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log 1*296000 6*11260 

log 

2*56260 

Jog T d 3 ,]t 9746 

log T ,f 

3-19746 

log nt" 2-91514 

T“ = iS 75 d-6 5 

m" - 822"-5 

ii. Calculation of Longitude of Perihelium and Mean Longi¬ 
tude. (48) to (54).— 

log tgT> 

9 ' 3 i 9 i 6 

O / ii 

= [92 5 8 

log cos i 9*98851 

= 171 8 29 

log tgu I 

9 ' 33 o 6 S 

363 13 37 

i+e = 1*24744 

— v z = 49 16 46 

1 — 

■e = 0-75256 

n = 52 30 23 

log tg ^ 

9*66i5i(w) log e 

9*39347 *0 = —0-684489 

log s/i — e 

9*93827 log sin w 

9-80091 (w) — (s 22 )= 0-156452 

c. logs/i+e 9-95199 log (s 22 ) 

9* 1 9438(m) M = —0-528037 

. w 
log tg — 

9 * 55 I 77 W ~' = 

—19 0 36' 33" *o= -39 0 13' 6" 
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E, = 1804, Oct. 5-458644 log M. 9-72263 (n) 

E 0 = 1805, Jan. o-oooooo log 206265" 5-31442 


T d = 86 d -54i356 log(s 23 


5 -° 37 ° 5 (m) 


log r 4 t'93722 
log to" 2-91514 

lo g (4 4 ) 4-85236 


(* 23 ) = -108905 
(*24)= 71180 

L" = - 37725 


L = —10 28 45 
11= 52 30 23 

L 0 = 41 58 22 


Syra, Greece : 1903 August 20. 


Ephemeris for Physical Observations of Saturn , 1903-4. 

By A, C. D. Crommelin. 

The discovery by Professor Barnard of a conspicuous white 
spot on Saturn last July has reawakened interest in the question 
of the planet’s rotation, and has illustrated the remarkable varia¬ 
tion of rate that prevails in different latitudes. At Mr. Den¬ 
ning’s suggestion I have computed an ephemeris giving the 
longitudes of central meridian and the times of transit of the zero 
meridian on two assumptions as to the rotation period, viz. : 

System I. io h I4 m *4 for equatorial spots (Hall). 

System II. io h 37 m '92 for the spot discovered by Barnard 
last July (period provisionally deduced by Denning). 

The quantities corresponding to P, B, B in the Jupiter 
ephemeris may be taken from the Nautical Almanac ephemeris of 
the ring, or more conveniently from the Connaissance des Temps , 
as this is fuller and contains the longitude of the Earth in the 
plane of the ring. In order to utilise the data of the Connais¬ 
sance without interpolation, I give the longitude of the central 
meridian for Paris midnight, corresponding to n h 5o m, 65 G.M.T.; 
but the transits of zero meridian are given for G.M.T. The 
quantities are to be interpolated for the time for which they are 
required, the equation of light having been already applied. The 
approximate increase of the longitude of central meridian in five 
days is 4219 0 in System I., 4063° in System II. I have com¬ 
menced. a new system of zero meridians, as Mr. Marth used a 
different period of rotation (io h i3 m *6) in his ephemerides, of 
which the last appeared in vol. lv. p. 164. The small figures 
between the transits of zero meridians denote the number of 
intervening rotations, and the table at the end facilitates the 

determination of the times of intervening transits. 
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